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Under  the  grant  N00014-89-J-1119  we  have  made  significant  progr^s 
in  a  number  of  areas  that  have  led  to  increased  understanding  of 
dynamics  in  solids  and  flames  and  in  the  development  of  new  experimental 
techniques  and  associated  theory  for  the  investigation  of  condensed 
matter  and  gas  phase  systems.  This  work  will  be  briefly  summarized  in 
this  section.  A  list  of  publications  that  have  resulted  from  the 
research  supported  by  this  grant  follows. 


A.  Surface  Selective  Transient  Grating  Experiments: 

Heat  Flow  at  Interfaces 

The  laser  induced  transient  grating  technique  which  has  been 
expanded  in  areas  of  application  and  sophistication  over  the  last  decade 
has  been  successfully  applied  to  the  investigation  of  a  wide  variety  of 
transient  processes  in  liquids,  solids,  and  gases.  Previous  theoretical 
treatments  of  transient  gratings  [1]  were  developed  by  direct  analogy 
with  the  diffraction  from  transparent  volume  holograms  [2,3]  where  the 
general  method  of  slowly  varying  amplitudes  was  utilized.  The  grating 
was  assumed  to  be  constant  in  amplitude  throughout  the  sample  and  the 
probe  beam  was  assumed  to  be  weakly  absorbed.  Recent  experimental 
investigations  on  samples  with  very  high  optical  densities  [4],  e.g. 
metals  or  absorbing  dielectrics  with  spatially  nonuniform  gratings, 
demonstrated  that  further  theoretical  development  is  required  in  which 
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the  effects  of  surfaces  are  included.  Recent  work  in  our  laboratory  has 
provided  such  a  formalism  for  describing  both  reflection  and 
transmission  geometries  of  diffraction  where  surface  effects  are 
specifically  considered  for  samples  with  or  without  a  high  optical 
density  and  an  arbitrary  grating  excitation  profile.  [5]  The  essential 
features  of  the  transient  grating  diffraction  efficiency  are  given  by 
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for  the  transmitted  geometry.  The  probe  wave-vector  in  vacuum  is 
written  and  the  grating  amplitude  or  peak  to  null  change  in  the 
dielectric  constant  is  denoted  A« .  The  reflected  diffracted  geometry 
(Eq.  1)  has  an  exponential  term  that  localizes  the  sampling  region  near 
the  incident  interface  while  the  transmitted  geometry  is  simply  the 
average  of  the  grating  amplitude  throughout  the  bulk. 

By  analyzing  different  diffraction  geometries,  the  transient 
grating  method  can  provide  information  on  the  spatial  distribution  of 
material  excitations.  This  is  illustrated  with  experimental  results 
from  a  thin  film  (220nm)  of  the  high  temperature  superconductor 
YBa2Cu20y_j^  (YBCO)  on  a  MgO(lOO)  substrate.  Absorption  of  the  561  nm 
laser  grating  excitation  beams  by  the  superconductor  produces  electronic 
excitations  that  rapidly  damp  (-Ips)  by  the  generation  of  a  large  number 
of  phonons.  This  incoherent  jump  in  the  phonon  population  is  heat,  and 
the  result  is  a  sinusoidal  heat  pattern.  The  sinusoidal  heating 
generates  a  sinusoidal  density  variation,  and  since  the  index  of 


refraction  is  density  dependent  (opto-elastic  constant),  a  holographic 
grating  is  produced. 
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The  heat  undergoes  diffusion  and  the  impulsive  heat  deposition 
drives  multiple  acoustic  waveguide  modes  [6].  The  probe  beam  is 
temporally  delayed  and  diffracted  from  the  transient  holographic 
grating.  Four  different  transient  grating  geometries  are  shown  in  Fig. 

1.  When  one  looks  at  the  reflected  geometries  (denoted  refl.  in  Fig. 

1) ,  the  probe  will  selectively  sample  a  region  near  the  incident  probe 
interface.  In  contrast,  the  transmission  geometries  (denoted  trans .  in 
Fig.  1)  will  sample  the  entire  sample  uniformly.  The  transient  grating 
method  can  therefore  directly  and  simultaneously  provide  information  on 
transient  processes  at  a  surface  or  interface  and  inside  the  bulk  of  a 
material . 

Figure  2  shows  the  results  of  the  transient  grating  experiments  at 
300  K.  The  solid  lines  thru  the  data  are  theoretical  fits  where  three 
dimensional  diffusion  is  assumed.  The  MgO  is  an  infinite  heat  sink 
since  the  thermal  conductivity  of  MgO  is  several  orders  of  magnitude 
higher  than  YBCO.  To  properly  fit  the  data,  a  barrier  limited  heat  flow 
model  must  be  used  at  the  MgO-YBCO  interface.  Several  different  models 
were  considered.  One  was  a  barrier  region  of  finite  thickness  inside 
the  YBCO  while  another  was  that  the  temperature  at  the  interface  was  a 
constant.  Neither  of  these  models  could  simultaneously  fit  all  four 
data  sets  without  individually  adjusting  the  thermal  constants  between 
geometries.  Using  a  surface  conductivity  model  where  the  temperature 
gradient  at  the  interface  Is  proportional  to  the  temperature  difference 
across  the  interface,  it  was  possible  to  simultaneously  fit  all  four 
geometries.  These  fits  are  the  ones  given  in  Fig.  2. 
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It  is  important  to  note  that  it  was  necessary  to  use  multiple 
geometries  to  properly  model  the  dynamics  inside  the  material,  at  the 
surface,  and  at  the  interface.  In  general,  thin  films  or  boundary 
layers  can  have  different  properties  than  those  of  the  bulk  material. 
Since  a  variety  of  important  physical  phenomena  occur  at  surfaces  and 
interfaces  and  many  new  devices  are  based  on  thin  film  technology,  it  is 
important  to  have  a  technique  that  can  nondestructively  measure  the 
properties  of  thin  films.  Conventional  techniques  for  measuring  thermal 
propagation  in  materials  rely  on  sensors  that  must  make  intimate  contact 
with  the  sample  without  disturbing  the  thermal  flow  itself  [7,8].  This 
requires  either  large  specially  shaped  samples  or  very  small  sensors. 

In  either  case  the  sample  size  is  limited  to  dimensions  much  larger  than 
100  nm.  Since  the  transient  grating  method  described  here  is  a  remote 
sensing  technique,  the  measurement  process  does  not  require  physical 
contact  with  the  sample.  Thus  the  sample  size  need  not  be  constrained 
to  dimensions  larger  than  those  of  the  detectors  and  very  thin  materials 
can  be  studied  down  to  atomic  dimensions.  This  illustrates  a  powerful 
new  way  to  study  thin  film,  interface,  and  surface  properties  where  bulk 
and  surface  effects  can  be  measured  separately.  While  the  results 
described  here  investigate  the  important  problem  of  interfacial  heat 
flow,  a  wide  variety  of  other  phenomena,  e.  g. ,  dynamics  of  molecules  or 
concentrations  of  chemical  species,  can  be  studied. 
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B.  Non- linear  Optical  Studies  of  Flames  and  Gas  Phase  Dynamics 

1.  Picosecond  Four-Wave  Mixing  Experiments  In  Sodium-Seeded  Flames 

The  ability  to  non- invasive ly  probe  the  local  properties  of  flames 
(e . g. , transport  properties,  species  concentrations,  temperature, 
reaction  rates,  and  collision  cross-sections)  over  small  (sub¬ 
millimeter)  distance  scales  is  essential  to  the  understanding  of 
combustion.  The  advent  of  lasers  has  allowed  some  of  these  processes  to 
be  studied,  and  techniques  such  as  Laser-Induced  Fluorescence  (LIF) , 
Raman  spectroscopy,  and  absorption  measurements  have  become  the 
mainstays  of  combustion  research  [9].  These  experiments  tend  to  be 
difficult  to  perform  however,  due  to  a  combination  of  weak  signals  and  a 
large  amount  of  background  noise  (both  from  the  flames  and  the  laser 
beams  themselves).  For  these  reasons,  four-wave  mixing  (FWM)  techniques 
have  proven  quite  useful  in  combustion  research;  these  techniques 
generally  provide  relatively  strong  signals  that  are  spatially  separated 
from  the  laser  beams.  Coherent  Anti-Stokes  Raman  Spectroscopy  (CARS) 
[10]  and  other  FWM  techniques  [11]  have  proven  useful  in  obtaining 
information  on  steady- state  flame  properties  such  as  concentration 
profiles.  Four-,  six-,  and  eight-wave  mixing  have  also  been  used  to 
study  collision-enhanced  spectra  in  sodium-seeded  flames  [12].  To  study 
the  fast  time-scale  dynamics  of  flames  (such  as  collisions,  velocity 
distributions,  and  reactions),  it  is  necessary  to  perform  experiments 
with  subnanosecond  laser  pulses . 

The  picosecond  transient  grating  [13]  is  a  unique  and  powerful 
tool  for  probing  fast  time-scale  dynamics.  Previous  work  in  this 
laboratory  has  demonstrated  that  the  transient  grating  is  an  effective 
method  for  probing  velocity  distributions  in  low-pressure  gases  [14]. 
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We  have  successfully  applied  the  transient  grating  to  the  study  of 
various  Na-seeded  pre-mixed  flames  [15]. 

Figure  3a  shows  an  intensity  grating  (a  grating  in  which  both 
excitation  beams  are  of  the  same  polarization)  decay  in  a  Na-seeded, 
methane/air  flame .  The  decay  is  exponential  with  a  decay  constant  of 
800  psec,  which  is  much  shorter  than  the  16  nsec  Na  excited-state 
lifetime  [16],  This  decay  is  insensitive  to  the  grating  fringe  spacing, 
and  therefore  corresponds  to  the  Na  excited- state  quenching  rate. 

Figure  3b  shows  a  polarization  grating  (a  grating  in  which  the 
excitation  beams  are  of  orthogonal  linear  polarizations)  decay  taken 
under  the  same  conditions  in  the  same  flame.  The  differences  between 
the  two  decays  is  striking.  Whereas  the  intensity  grating  decay  is 
smooth,  the  polarization  grating  decay  shows  large  oscillations  at  the 
Na  ground-state  hyperfine-splitting  frequency.  More  remarkably,  the 
decay  time  of  the  polarization  grating  is  seven  times  longer  than  that 
of  the  intensity  grating  at  this  fringe  spacing.  Moreover,  the 

polarization  grating  decay  constant  is  sensitive  to  fringe  spacing. 

2 

Figure  4  is  a  plot  of  this  constant  vs.  (Zw/d)  ,  where  d  is  the  grating 
fringe  spacing.  The  points  in  this  plot  fall  along  a  straight  line 
whose  slope  is  related  to  the  Na  diffusion  constant.  Using  these 
techniques ,  we  have  been  able  to  measure  four  picosecond  d3mamic 
processes  in  various  Na-seeded  flames:  excited-state  quenching,  excited- 
state  population  scattering,  ground-state  magnetic-sublevel  population 
scattering,  and  the  diffusion  of  Na  atoms.  We  will  soon  be  using  these 
techniques  to  probe  the  dynamics  of  intrinsic  flame  species. 
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2.  The  Grating  Decomposition  Method 

Four-wave  mixing  (FWM)  techniques  have  become  increasingly 
powerful  tools  for  probing  phenomena  in  all  phases  of  matter.  An 
important  special  case  of  FWM,  the  transient  grating  [17]  (TG) ,  has 
proven  useful  in  studying  a  multitude  of  interesting  processes. 

Although  the  transient  grating  is  a  versatile  technique,  TG  data  are  not 
always  simple  to  interpret,  especially  those  arising  from  orientational 
gratings  (gratings  in  which  the  direction  of  the  electric  field  changes 
along  each  fringe;  an  example  of  this  is  the  linear  polarization  grating 
[18],  in  which  the  excitation  beams  are  cross-polarized).  Historically, 
two  approaches  have  been  taken  in  the  interpretation  of  orientational  TG 
data.  The  first  approach  makes  predictions  based  on  the  spatial 
dependence  of  the  amplitude  and  direction  of  the  grating  electric  field 
[19,20].  While  this  approach  has  the  advantage  of  being  simple  and 
intuitive,  it  provides  little  calculational  power.  Even  worse,  electric 
field  pictures  may  at  times  lead  to  unphysical  predictions,  such  as 
higher-order  diffraction  at  low  excitation  beam  intensities. 

The  other  approach  that  has  been  used  to  understand  TG  experiments 
and  FWM  in  general  [21-23]  is  the  density  matrix  formalism,  usually  in 
the  form  of  diagrammatic  perturbation  theory.  In  this  formalism, 

Feynman  diagrams  are  used  to  depict  integrals  that  describe  the 
couplings  of  material  states  by  laser  fields  and  the  evolution  of  these 
material  states  in  between  field  interventions.  A  weakness  of  this 
theory  as  it  is  usually  applied  to  orientational  TG  experiments, 
however,  is  that  spatial  information  about  the  excitation  processes  is 
not  available  in  a  form  that  provides  physical  insight.  This 
information  can  be  essential  to  the  understanding  and  calculation  of 
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orientational  gratings  and  their  decays  for  several  reasons.  First, 
decay  constants  must  be  added  in  a  completely  ad  hoc  manner;  in  some 
systems  (such  as  gas-phase  sodium  atoms),  the  deca>  constants  for  both 
linear  intensity  gratings  and  linear  polarization  gratings  are 
calculationally  identical,  even  though  the  experimentally  observed  decay 
constants  may  differ  greatly.  Such  differences  arise  because  of  the 
different  spatial  dependences  of  the  excitation  processes  in  each 
experiment.  Second,  this  theory  cannot  treat  secondary  effects  such  as 
thermal  gratings  [24]  or  molecular  alignment  (such  as  through  the 
nuclear  Optical  Kerr  Effect  [20]),  even  though  these  effects  may  be 
major  contributors  to  grating  signals.  It  is  difficult  to  extend 
diagrammatic  perturbation  theory  to  treat  such  processes  without 
understanding  the  spatial  nature  of  the  excitation  process. 

We  have  shown  that  it  is  possible,  through  the  decomposition  of 
excitation  beam  polarizations  into  various  component  polarizations,  to 
transform  diagrammatic  perturbation  calculations  of  orientational 
gratings  into  sums  of  linear  (±45°)  and  circular  (right  and  left) 
intensity  grating  calculations  [25]  (see  Fig.  5).  In  each  of  these 
component  int;ensity  gratings,  the  two  excitation  beams  have  the  same 
polarization,  and  therefore  spatial  information  is  restored.  Thus 
spatial  information  (and  therefore  insight)  is  included  in  orientational 
grating  calculations.  We  have  also  introduced  effective  two- interaction 
matrix  elements  (ETIMEs) ,  which  simplify  some  of  the  summations  that 
must  be  performed  in  diagrammatic  perturbation  calculations.  We  have 
used  the  GDM  and  ETIMES  to  simplify  the  treatment  of  picosecond 
absorption  and  emission  orientational  gratings  and  nuclear  Optical  Kerr 
Effect  orientational  gratings  in  isotropic  molecular  liquids.  The  GDM 
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approach  to  the  latter  problem  clears  up  inconsistencies  that  arose  in 
previous  treatments.  In  addition,  we  have  successfully  used  the  GDM  and 
ETIMEs  to  treat  a  problem  that  cannot  be  properly  treated  within  the 
confines  of  standard  diagrammatic  perturbation  theory:  orientational 
gratings  of  gas-phase  Na  atoms.  We  have  been  able  to  accurately  model 
the  behavior  of  polarization  gratings  on  the  Na  D  lines  in  low-pressure 
gas  cells  and  in  Na- seeded  flames  [26] ,  while  at  the  same  time  greatly 
reducing  the  necessary  calculations.  The  GDM  and  ETIMEs  will  find  many 
applications  in  orientational  grating  problems. 

3. Picosecond  Time-Scale  Phase-Related  Optical  Pulses 

The  development  of  optical  coherence  techniques  for  probing  the 
properties  of  systems  of  atoms  and  molecules  has  in  many  instances 
closely  followed  the  path  taken  by  nuclear  magnetic  resonance  (NMR) 
techniques.  In  general,  however,  optical  pulse  sequences  that  parallel 
magnetic  resonance  have  not  been  able  to  employ  one  important  property 
of  the  pulses  in  most  NMR  sequences:  well-defined  phase  relationships 
[27,28].  Without  the  ability  to  control  the  relative  phases  of 
ultrafast  (picosecond  or  femtosecond)  laser  pulses  that  are  separated  in 
time,  it  is  impossible  to  utilize  many  potentially  useful  experimental 
techniques,  such  as  the  fluorescence-detected  photon  echo  [29]  and  the 
dipolar  echo  [30]. 

We  have  developed  a  technique  that  can  potentially  be  used  to 
phase  lock  any  arbitrary  number  of  ultrafast  pulses  [31].  In  addition, 
this  method  can  sef/e  as  an  add-on  to  existing  ultrafast  laser  systems. 
We  have  used  this  technique  to  phase -lock  a  pair  of  picosecond  time 
scale  laser  pulses  in  order  to  measure  the  optical  coherence  decay  of 
the  line  of  sodium  vapor. 
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In  this  experiment,  a  40-psec,  phase -locked,  temporally  separated 
laser  pulses  that  are  resonant  with  the  Na  transition  are  passed 
through  a  Na  gas  cell .  The  resultant  incoherent  fluorescence  is 
monitored  at  a  right  angle  to  the  path  of  the  laser  beams.  As  the 
relative  phase  between  the  pulses  is  scanned,  the  fluorescence  from  the 
Na  cell  is  sinusoidally  modulated.  Such  a  scan  is  illustrated  in  Fig. 

6,  in  which  the  delay  time  is  107  psec.  By  plotting  the  peak-null 
difference  in  the  fluorescence  oscillations  as  a  function  of  the  delay 
time  between  the  pulses,  we  were  able  to  plot  the  coherence  decay  of  the 
Na.  This  data  is  illustrated  in  Fig.  7,  along  with  a  theoretical  fit 
that  uses  no  adjustable  parameters.  The  fluorescence  oscillations  in 
Fig.  6  and  the  excellent  fit  between  theory  and  experiment  in  Fig.  7 
illustrate  that  the  relative  phase  between  the  laser  pulses  is  indeed 
well-defined.  This  method  of  phase -locking  promises  to  be  useful  in 
ultrafast  experiments  in  all  phases  of  matter. 

C.  Ultrafast  Dynamics  in  Complex  Liquids 

The  rotational  dynamics  of  liquids  have  been  under  intense 
investigation  in  recent  years.  There  exists  a  large  number  of 
techniques  for  studying  these  dynamics,  including  fluorescence 
depolarization,  dynamic  light  scattering  ,  optical  Kerr  effect  (OKE) , 
nuclear  magnetic  resonance,  electron  spin  resonance,  and  Raman 
scattering.  The  majority  of  the  time  domain  experiments  have  been  done 
on  very  small,  simple  molecules  such  as  CS2  [32,33].  These  experiments 
were  primarily  done  to  investigate  the  ultrafast  librational  dynamics. 
The  rotational  diffusion  of  these  molecules  has  usually  been  described 
by  a  single  exponential. 
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The  rotational  dynamics  of  molecules  has  been  assumed  to  be 
adequately  described  by  a  modified  Debye -Stokes -Einstein  (DSE)  theory. 
This  theory,  which  was  developed  to  describe  the  d3mamics  of  a  rotating 
sphere  in  a  featureless  continuum,  relates  the  rotational  diffusion  time 
Trot  bulk  viscosity  of  the  medium,  r) ,  and  the  temperature,  T 

[34]. 

’’rot  “ 

The  proportionality  constant,  is  the  effective  volume  of  the 

molecule.  Molecular  reorientational  motions  that  are  described  by  DSE 
theory  can  be  referred  to  as  hydrodynamic. 

The  ultrafast  component  of  the  rotational  dynamics  of  liquids  can 
not  be  described  by  DSE  model.  At  times  shorter  than  1-2  ps ,  the 
dynamics  are  dominated  by  librations  and  collision- induced  (Cl) 
phenomena.  The  equilibrium  ensemble  of  librations  is  altered  by  the 
fsec  laser  pulse.  The  broad  bandwidth  of  the  fsec  laser  pulse  permits 
coherent  Raman  excitation  of  librations  of  frequencies  within  the  pulse 
bandwidth  [35]. 

The  experiments  that  are  described  here  were  performed  with  a 
subpicosecond  laser  system  developed  in  our  labs  [36] .  The  output  of  a 
cw  mode -locked  Nd;YAG  laser  was  pulse -compressed  and  frequency  doubled 
to  s3mchronously  pump  a  dye  laser,  providing  tunable  300  fs  pulses  at  82 
MHz.  A  Q-switched,  and  mode-locked  Nd:YAG  laser  is  cavity  dumped  at 
1.75  kHz  and  then  frequency  doubled.  The  700  fiJ  single  pulses  at  532  nm 
are  used  to  amplify  the  dye  laser  output.  Transform- limited  150-300  fs 
amplified  pulses  of  energies  up  to  50  fiJ  were  obtained  in  the  wavelength 
range  of  560-680  nm  (R6G  and  DCM  dyes). 

The  OKE  was  studied  using  a  transient  grating  technique,  first 
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used  by  Eyring  and  Payer  [37].  The  transient  grating  approach  has  a 
number  of  advantages  over  a  pump-probe  transient  birefringence  method, 
including  background- free  detection  and  the  ability  to  separate  the 
electronic  and  nuclear  contributions  to  the  OKE.  As  first  described  by 
Etchepare  et  al  [38],  and  extended  by  Deeg  and  Payer  [39],  physical 
processes  having  different  x(3)  symmetries  can  be  observed  independently 
by  choosing  the  correct  polarizations  for  the  three  input  beams  and  the 
signal.  The  nuclear  contribution  can  be  eliminated,  providing 
observation  of  the  electronic  response  (which  gives  the  instrument 
response).  If  the  electronic  contribution  is  eliminated,  direct 
observation  of  the  nuclear  response  is  permitted. 

The  reorientational  relaxation  of  biphenyl  showed  relatively 
straight  foirward  behavior  [40].  At  times  greater  than  2  ps ,  the 
reorientational  dynamics  of  biphenyl  can  be  described  by  a  bi¬ 
exponential  decay.  The  slower  component  corresponds  to  the  tumbling 
motion  of  the  biphenyl  molecule,  whereas  the  faster  component  is 
associated  with  the  torsional  motion  of  the  two  phfenyl  rings  around  the 
central  C-C  bond.  A  temperature  dependence  showed  hydrodynamic  behavior 
for  the  66  °C  temperature  range  studied.  The  ultrafast  dynamics  (<  2 
ps)  was  described  as  a  pair  of  overdamped  oscillators,  one  corresponding 
to  each  of  the  rotational  diffusion  decays.  Explicit  linkage  of  the 
librational  and  rotational  diffusion  decays  was  made  for  the  first  time. 
Typical  data  are  shown  in  Figure  8. 

A  similar  study  of  the  isotropic  phase  of  the  nematic  liquid 
crystal  pentylcyanobiphenyl  (5CB)  showed  vastly  more  complex  behavior 
[41-43] .  There  is  a  very  slow  component  to  the  decay  (up  to  100  ns  near 
the  phase  transition)  corresponding  to  the  motion  of  the  microscopic 


TG— Signal  intensity  (arb.  units)  TG— Signal  intensity  (arb.  units) 


-2  0  2  4 

time  (ps) 


Figure  8 


0 


13 


psuedonematic  domains .  The  temperature  dependence  of  this  component  of 
the  decay  is  well  described  by  Landau-deGennes  theory  [13] 

Tj.  a  f?/(T-T*)  (4) 

where  T  is  the  virtual  phase  transition  temperature.  This  model  fits 
the  data  up  to  30-35  °C  above  the  phase  transition.  At  this 
temperature,  the  size  of  the  psuedonematic  domains  becomes  of  the  order 
of  a  few  molecules.  In  that  same  temperature  range,  all  of  the  faster 
dynamics  are  temperature  independent.  This  is  a  breakdown  of  DSE  theory 
and  represents  a  decoupling  of  the  d3mamics  from  the  bulk  viscosity.  As 
can  be  seen  in  fig.  9,  a  temperature  dependence  in  the  slowest  non¬ 
collective  component  develops  at  higher  temperatures  (i.e.,  when  the 
psuedonematic  domains  become  vanishingly  small) .  When  5CB  is  put  in  n- 
heptane  solutions,  the  collective  component  of  the  decay  disappears,  and 
the  remaining  components  behave  hydrodynamically . 

A  study  of  2 -ethyl  naphthalene  shows  complex  behavior  in  a  simple 
liquid  [45].  A  times  greater  than  2  ps  (after  the  ultrafast  librational 
dynamics  have  damped  out) ,  the  decay  can  be  accurately  described  by  a 
triple  exponential  with  each  component  having  4  to  6  factors  of  e  decay. 
Over  the  80  °C  temperature  range  studied,  the  slowest  two  components 
exhibited  hydrodynamic  behavior.  The  fastest  component  was  temperature 
independent  from  2-40  °C.  In  analogy  to  the  collective  motions  in  5CB, 
this  represents  the  domination  of  the  fast  dynamics  by  local  structures 
in  the  liquid.  These  local  structures  have  dynamics  that  depend  on  the 
details  of  the  local  intermolecular  interactions  rather  than  on  the  bulk 
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D.  Probing  Disordered  Solids  with  Variable  Time  Scale  Optical  Dephaslng 
Measurements 

Dynamics  and  interactions  in  glasses  and  other  disordered 
materials  such  as  imperfect  crystals  or  proteins  can  be  studied  with 
optical  line  narrowing  techniques  applied  to  chromophores  embedded  in 
the  disordered  material.  Prior  to  the  use  of  optical  methods,  a  wide 
variety  of  classical  methods,  such  as  heat  capacities,  were  applied  to 
the  study  of  glasses.  Results  have  been  successfully  interpreted  in 
terms  of  a  model  of  the  glass  potential  surface  described  in  terms  of 
the  Two  Level  Systems  (TLS)  [46,47]  model.  Because  of  the  very  large 
extent  of  inhomogeneous  broadening  in  glassy  systems,  it  is  necessary  to 
apply  line  narrowing  methods  to  provide  information  on  the  dynamics  and 
interactions  of  atoms  and  molecules  with  their  environments. 

A  number  of  optical  line  narrowing  methods  have  been  developed  to 
remove  inhomogeneous  broadening.  Among  these  are  hole  burning, 
fluorescence  line -narrowing,  accumulated  grating  echoes,  stimulated 
photon  echoes,  and  photon  echoes.  The  photon  echo  and  the  stimulated 
photon  echo  are  the  direct  optical  analogs  of  the  magnetic  resonance 
spin  echo  and  stimulated  spin  echo.  In  magnetic  resonance  it  is  known 
that  the  spin  echo  measures  the  homogeneous  dephasing  time,  T2,  and  that 
the  stimulated  echo  can  have  additional  contributions  to  the  dephasing 
time  from  spectral  diffusion  occurring  on  time  scales  slower  than  T2. 
Until  recently  hole  burning  and  fluorescence  line  narrowing  have  been 
the  line  narrowing  techniques  applied  to  the  study  of  glasses.  Various 
workers  have  defined  the  observables  of  these  techniques  as  the 
homogeneous  linewidth.  Each  technique,  however,  is  sensitive  to 
processes  over  a  different  range  of  time  scales  and,  in  general,  will 
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yield  different  linewidths  [48-50]. 

The  line  narrowing  experiments  have  been  treated  as  if  they  are 
absorption  experiments  using  optical  absorption  formalisms.  The  line 
narrowed  lineshape  was  taken  to  be  related  to  the  Fourier  transform  of 
the  two-time  transition  dipole-moment  correlation  function.  However, 
the  two  time  correlation  function  is  insufficient  to  describe  actual 
line -narrowing  experiments,  such  as  hole-burning,  that  are  performed  on 
glasses  and  crystals  [48].  In  fact  the  two  time  correlation  function 
description  of  experiments  has  been  misleading.  It  indicated  that  all 
of  the  line  narrowing  experiments  measured  the  same  observable,  the 
homogeneous  line.  It  caused  the  fundamentally  distinct  nature  of 
different  experiments,  i.  e.,  the  characteristic  time  scale  associated 
with  each  experiment,  to  be  absent  from  the  theoretical  treatments. 

Using  the  correct  description  of  the  different  experiments  makes  it 
possible  to  obtain  information  on  the  broad  distribution  of  rates  of 
dynamical  processes  in  complex  systems  [48-50] . 

A  theoretical  formalism  has  been  developed  by  Berg  et  al.  [51]  and 
extended  by  Bai  and  Payer  [52]  that  relates  the  experimental  observables 
in  the  various  types  of  line  narrowing  experiments  to  details  of  the 
djmamics  and  intermolecular  interactions  in  complex  systems  such  as 
glasses.  Each  line -narrowing  technique  can  be  regarded  as  a  specific 
case  of  a  four-wave  mixing  experiment  in  which  three  input  beams  induce 
and  subsequently  interfere  with  a  polarization  in  a  medium  to  generate  a 
fourth  (signal)  beam.  The  frequency  domain  experiments  have  specific 
time  domain  analogs  obtained  through  Fourier  transforms.  The  differences 
between  the  techniques  are  reduced  to  the  timing  between  the  pulses  in 
the  time  domain  description.  It  has  also  been  shown  that  all  of  the 
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line -narrowing  experiments  currently  being  performed  on  glasses  are 
describable  in  terms  of  correlation  functions  that  are  variants  of  the 
stimulated  echo  four  time  correlation  function.  In  a  stimulated  echo, 
three  pulses  are  applied  at  times,  0,  r  and  r+T^.  A  fourth  pulse,  the 
signal,  emerges  at  time  2r+T^.  The  correlation  function  is  usually 
specified  by  the  time  intervals  between  the  pulses,  i.e.,  r,  T^,  and  r. 
The  detailed  analysis  shows  that  the  differences  among  the  various  line 
narrowing  experiments  is  the  duration  of  the  interval  T^,  which  has  a 
major  influence  on  the  observables  in  the  experiments  [48-52]. 

First  consider  the  photon  echo  experiment.  It  is  the  limit  of  a 
stimulated  echo  experiment  with  T^-0.  The  second  and  third  pulses  in 
the  sequence  are  time  coincident.  Inhomogeneous  broadening  is 
eliminated  by  the  echo  pulse  sequence.  The  polarization  in  the  sample 
decays  as  the  separation  between  the  two  excitation  pulses  is  increased. 
If  the  homogeneous  line  shape  is  Lorentzian,  the  polarization  decays  as 
exp[-2r/T2]  and  the  signal  intensity  decays  as  exp[-4r/T2]  where  T2  is 
the  homogeneous  or  transverse  dephasing  time.  In  addition  to  the 
dephasing  induced  by  dynamics,  population  decay  also  contributes  to  the 
coherence  loss.  The  pure  dephasing  time,  T2  can  be  obtained  by 
removing  the  population  decay  contribution. 

In  a  stimulated  echo,  a  pulse  excites  the  chromophores  as  before 
but  at  time  r ,  the  second  pulse  is  split  into  two  pulses  separated  by 
the  interval  T^.  The  echo  appears  as  a  fourth  pulse.  The  stimulated 
photon  echo  is  sensitive  to  homogeneous  dephasing  during  the  interval  r 
between  the  first  and  second  pulses  and  between  the  third  and  echo 
pulses.  It  is  also  sensitive  to  frequency  fluctuations  (spectral 
diffusion)  that  occur  on  the  slower  time  scale,  T^. 


17 


Hole  burning  is  a  frequency  domain  line  narrowing  experiment.  In 
a  hole  burning  experiment,  a  narrow  band  laser  excites  a  narrow 
transition  frequency  subset  of  the  chromophores  in  the  inhomogeneous 
line.  A  spectrum  bracketing  the  point  of  excitation  will  show  reduced 
absorption,  the  hole.  Molecules  that  were  excited  by  the  laser  can 
contribute  to  the  hole,  these  being  ones  either  initially  on  resonance 
with  the  laser,  or  ones  that  had  fluctuations  in  energy  during  the  hole 
burning  period  that  brought  them  into  resonance  with  the  laser.  In  the 
time  interval  between  burning  and  reading  the  hole,  fluctuations  in 
transition  energies  can  bring  molecules  into  the  frequency  region  of  the 
hole.  Therefore,  fluctuations  in  frequency  that  occur  on  the  time  scale 
from  burning  through  reading  the  hole  will  determine  the  hole  width. 
Since  this  time  scale  is  typically  long  (seconds),  the  experiment  will 
be  sensitive  to  very  slow  dynamical  processes  (very  slow  spectral 
diffusion).  It  has  been  formally  proven  that  hole  burning  is  the 
frequency  domain  equivalent  of  the  stimulated  echo  experiment  [51,52]. 
The  spectrum  of  the  hole  is  determined  by  the  Fourier  transform  of  the 
stimulated  echo  correlation  function.  Hole  burning  is  a  frequency 
domain  stimulated  echo  experiment.  In  a  hole  burning  experiment,  T^  is 
essentially  the  time  required  uo  do  the  experiment,  from  burning  through 
reading.  The  typical  time  scale  for  a  photon  echo  experiment  is  100  ps . 
Fluctuations  occurring  on  a  ^isec  or  msec  time  scale  are  static  on  the 
time  scale  of  a  photon  echo  experiment,  and  do  not  contribute  to  the 
measured  dephasing  time.  In  contrast,  fisec  and  msec  time  scale 
fluctuations  are  rapid  on  the  seconds  time  scale  of  a  hole  burning 
experiment,  and  contribute  to  the  measured  dephasing  time.  Therefore, 
echo  and  hole  burning  experiments  should  not  measure  the  same  dephasing 
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times  in  systems  which  have  broad  distributions  of  fluctuation  rates. 
Changing  T^,  the  time  scale  of  experiments,  provides  a  path  for  the 
examination  of  dynamics  over  a  very  wide  range  of  times  from  psec  to 
10,000  sec,  or  longer. 

The  dephasing  time  measured  by  photon  echoes  is  substantially 
longer  than  the  dephasing  time  measured  by  hole  burning  experiments 
[48] .  The  additional  broadening  of  a  hole  arises  from  slow  time  scale 
spectral  diffusion.  At  low  temperature,  the  ratio  of  the  hole  line 
width  to  echo  determined  line  width  is  proportional  to  the  integral  of 
the  fluctuation  rate  distribution,  P(R) ,  from  time  scale  of  the 
homogeneous  dephasing  time  to  the  time,  T^,  associated  with  the  hole 
burning  experiment  [52].  P(R)  is  the  probability  of  having  a  dynamical 
process  of  rate  R  in  the  glass.  Thus  if  T^  (time  between  burning  and 
reading  the  hole)  is  varied  over  a  range  of  times  in  which  some  rates 
are  active,  i.  e.,  1/R  falls  within  the  experimental  time  range,  the 
holewidth  will  change.  The  change  in  holewidth  can  be  used  to 
quantitatively  determine  P(R) . 

Recently  we  performed  experiments  of  this  nature .  We  examined 
cresyl  violet  in  ethanol,  ethanol-d,  and  glycerol  glasses  and  resorufin 
in  ethanol  [48-50].  For  example,  holes  were  measured  at  different  times 
after  burning  in  the  cresyl  violet/ethanol  system.  As  the  time  between 
burning  and  reading  increased,  the  hole  broadens.  For  cresyl  violet  in 
ethanol  at  1.30  K  the  holes  broadened  from  2.2  GHz  to  3.2  GHz  between 
100  msec  and  7000  sec.  The  holewidth  was  observed  to  vary  over  the 
entire  range  of  study,  but  the  change  in  width  is  most  pronounced  in  the 
range  of  50  seconds . 

By  performing  a  detailed  theoretical  analysis  of  the  four  time 
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correlation  function  for  hole  burning,  the  dependence  of  the  hole  width 
on  the  waiting  time,  T^,  was  obtained.  It  is  possible  to  determine  the 
fluctuation  rate  distribution  by  examining  the  holewidth,  or  the 
dephasing  time  in  time  domain  experiments,  as  a  function  of  T^.  A  log 
normal  distribution  (Gaussian  on  a  log  time  scale)  was  found  to  be  the 
form  of  P(R)  that  fit  the  data.  At  low  temperatures,  TLS  relaxation  is 
caused  by  phonon  assisted  tunneling.  For  this  case,  a  log  normal 
distribution  of  fluctuation  rates  is  equivalent  to  a  Gaussian 
distribution  of  tunneling  parameters. 

The  experiments  were  repeated  on  resorufin  in  ethanol,  and 
identical  results  were  obtained.  This  proves  that  the  d3mamics  in  this 
time  range  are  an  intrinsic  property  of  the  glass,  independent  of  the 
chromophore.  Experiments  on  glycerol  glass  covering  the  same  range  of 
T^  displayed  no  hole  broadening.  Again  demonstrating  that  the  dynamics 
observed  in  ethanol  glass  is  an  intrinsic  property  of  the  glass  and  not 
a  function  of  the  experimental  procedures. 

The  experiments  described  above  used  a  conventional  hole  burning 
procedure,  i.  e.,  the  hole  was  burned  and  following  a  time  T^,  the  hole 
was  read  out  by  sweeping  the  laser  frequency.  We  have  extended  the 
measurements  to  T^s  (time  between  burning  and  reading)  as  fast  as  10 
psec  [49] .  To  accomplish  this  we  measure  the  entire  time  response  at  a 
given  frequency  and  then  repeat  this  for  a  number  of  frequencies.  This 
avoids  the  necessity  of  sweeping  the  reading  frequency  across  the  hole. 
To  do  this  we  employ  a  300  MHz  acousto-optic  modulator  (AOM) .  The 
narrow  band  dye  laser  beam  is  split,  one  part  becomes  the  burning  beam. 
This  is  passed  through  two  AOMs  which  are  used  to  turn  it  on  and  off 
rapidly  to  make  the  burning  pulse.  The  other  part  of  the  dye  laser  beam 
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passes  through  the  300  MHz  AOM.  The  diffracted  beam  is  shifted  from  the 
burning  beam  by  exactly  300  MHz .  To  read  at  600  MHz ,  the  diffracted 
beam  is  passed  back  through  the  AOM.  The  beam  is  again  diffracted,  and 
is  now  shifted  600  MHz  from  the  burn  frequency.  This  process  can  be 
repeated  up  to  six  time  with  a  single  AOM  with  more  than  sufficient 
intensity  to  read  the  hole.  Thus  reading  frequencies  of  0,  300,  600, 
900,  1200,  1500,  and  1800  MHz  shift  from  the  center  frequency  are 
obtained.  This  can  be  done  on  either  side  of  the  line  by  taking  either 
the  +1  or  -1  diffraction  order. 

At  each  probe  frequency  the  full  time  dependence  is  measured.  The 
holes  used  are  transient,  presumably  triplet  bottle  neck  holes,  which 
refill  on  the  100  msec  time  scale  for  creysl  violet  in  ethanol. 
Therefore,  about  1000  shots  can  be  averaged  to  obtain  good  quality  data. 
Great  care  is  taken  to  test  for  saturation,  heating  and  other  unwanted 
side  effects. 

The  holes  at  10  /xsec  are  1.2  GHz  compared  to  3.2  GHz  measured  at 
7000sec.  Between  10  /isec  and  100  msec  the  increase  in  width  plotted  on 
a  log  time  scale  is  well  fit  by  a  straight  line,  indicating  a  P(R)  -  1/R 
distribution.  Therefore,  using  hole  burning  we  have  examined  hole 
widths  from  10  /isec  to  10,000  sec,  i.  e.,  over  9  decades  of  time. 

On  fast  time  scales  we  have  employed  photon  echo  and  stimulated 
echo  experiments  [48,50].  The  photon  echo  experiments  demonstrates  that 
the  fluctuation  rate  distribution  from  1  psec  to  1  nsec  also  has  the 
form  P(R)  -  1/R.  Extrapolation  of  the  msec  to  /isec  hole  burning  data  to 
the  nsec  time  scale  produces  essentially  perfect  agreement  with  the 
photon  echo  data.  This  indicates  that  a  1/R  distribution  exists  in 
ethanol  from  psecs  to  100  msec,  i.  e.,  over  11  decades  of  time. 
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As  discussed  above,  the  hole  burning  experiments  are  the  Fourier 
transform  of  the  stimulated  echo  experiment.  To  work  on  fast  time 
scales  it  is  preferable  to  work  in  the  time  domain  rather  than  the 
frequency  domain.  To  extend  the  short  time  echo  experiments  to  longer 
times,  to  avoid  the  necessity  of  doing  very  fast  hole  burning 
experiments ,  and  to  fill  in  the  gap  between  psec  time  scale  hole  burning 
experiments  and  psec  time  scale  photon  echo  experiments,  we  are  using 
the  stimulated  photon  echo.  Because  of  spectral  diffusion  on  a  fast 
time  scale,  the  stimulated  echo  decay  is  faster  than  the  two  pulse  echo 
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decay.  This  is  equivalent  to  the  classic  NMR  experiments  that  compare 
spin  echoes  and  stimulated  echoes  to  measure  spectral  diffusion.  The 
results  confirm  that  P(R)  is,  in  fact,  1/R  for  fast  fluctuation  rates. 

The  net  result  is  the  observation  of  hole  broadening  over  9 
decades  of  time  and  the  determination  of  the  low  temperature  ethanol 
glass  fluctuation  rate  distribution  over  16  decades  of  time.  The 
increase  in  observed  dynamical  line  widths  by  a  factor  of  approximately 
8  from  the  psec  time  scale  to  the  10,000  sec  time  scale  is  caused  by 
spectral  diffusion  induced  by  the  fluctuating  Two  Level  System  structure 
of  glasses.  The  data  is  analyzed  using  a  detailed  theoretical 
description  of  the  four  time  correlation  functions  that  describe  the 
various  line  narrowing  experiments. 
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E.  Theory  of  Shock  Induced  Multiphonon  Up-pumplng 

We  have  developed  [53]  a  theory  for  vibrational  up-pumping,  defect 
hot  spot  formation,  and  the  onset  of  chemistry  induced  by  a  shock  wave 
in  a  solid  composed  of  large  organic  molecules,  e.  g. ,  RDX.  The  model 
and  calculations  describe  the  flow  of  energy  in  a  shocked  solid.  The 
shock  excites  the  bulk  phonons,  which  rapidly  achieve  a  state  of  phonon 
equilibrium  characterized  by  a  phonon  quasi- temperature .  The  excess 
energy  subsequently  flows  into  the  molecular  vibrations ,  which  are 
characterized  by  a  vibrational  quasi- temperature .  The  multiphonon  up 
pumping  process  occurs  because  of  enharmonic  coupling  terms  in  the  solid 
state  potential  surface.  Of  central  importance  are  the  lowest  energy 
molecular  vibration,  or  "doorway"  modes,  through  which  mechanical  energy 
enters  and  leaves  the  molecules.  Using  reasonable  experimental 
parameters,  it  is  found  that  the  quasi- temperature  increase  of  the 
internal  molecular  vibrations  and  equilibration  between  the  phonons  and 
vibrations  is  achieved  on  the  time  scale  of  a  few  tens  of  picoseconds. 

A  new  mechanism  is  presented  for  the  generation  of  "hot  spots"  at 
defects.  Defects  are  postulated  to  increase  the  enharmonic  coupling  in 
a  localized  area  of  the  solid,  causing  the  vibrational  temperature  in 
this  area  to  briefly  overshoot  the  bulk.  The  influence  of  the  higher 
defect  vibrational  temperature  on  chemical  reactivity  is  estimated.  It 
is  shown  that  even  small  increases  in  defect  anharmonic  coupling  have 
profound  effects  on  the  probability  of  shock  induced  chemistry.  The 
anharmonic  defect  model  predicts  a  size  effect.  The  defect  enhanced 
chemical  reaction  probability  is  reduced  as  the  particle  size  is 


reduced. 
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